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Abstract 

Photoproduction of J/ ip is considered including the color-octet contributions 
from the various partial wave states, ^'^"'"^Lj = ^5o, ^Si and ^Pj. The produc- 
tion cross section depends on three new nonperturbative parameters defined 
in NRQCD, called the color-octet matrix elements. Using the color-octet 
matrix elements determined by fitting the J/ip production at the Tevatron, 
we find that the color-octet {cc)%{}So and ^Pj) contributions to the J/ijj 
photoproductions at the fixed target experiments and HERA are too large 
compared to the data on (7(7 + p ^ J/i' + X) in the forward direction, the 
z distribution of J/V'- The Pj, distribution of J/ijj and the total inelastic 
J /ijj production rate as a function of ^/s^ are predicted including color-octet 
contributions. We also briefly digress on the B — > J/ijj -\- X and observe 
the similar situation. This may be an indication that the color-octet matrix 
elements determined from the J/ip production at the Tevatron, especially 
(OlOg (^S'o)|0) and (OlCg (^Pj)|0), might have been overestimated by an order 
of magnitude. 
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I. INTRODUCTION 



In the conventional approach, the inelastic (inclusive) J /if) photoproduction has been 
studied in the framework of perturbative QCD (PQCD) and the color singlet model In 
this model, one considers '-f + g ^ J/ip + g which can produce high J/ip^s at the ep or 
7P collision. However, the same approach, when applied to the J/ip oi ip production at the 
Tevatron, severely underestimates the productions rate . In order to reconcile the data and 
PQCD predictions, a new mechanism for heavy quarkonium productions has been suggested 
0, the color-octet gluon fragmentation into J /ip. Also, the color-octet mechanism in heavy 
quarkonium productions at hadron colliders through the color-octet {cc)^ pair in various 
partial wave states "^^^^Lj has been considered beyond the color-octet gluon fragmentation 
approach [^. The main motivation is that inclusive T productions at the Tevatron also 
show the excess of the data over theoretical estimates of the productions based on PQCD 
and the color singlet model |^. Here, the px of the T is not that high so that the gluon 
fragmentation picture may not be a good approximation any more. In Refs. and 0, 
a large class of color-octet diagrams has been considered which can contribute to the J/ifj 
production at hadron colliders. At the partonic level, there appear new 2 — > 1 subprocesses 



gg-.(cc)e^f)), (1.1) 
gg^{cc){'S^^Ki'pf^), (1.2) 

at the short distance scale, and the subsequent evolution of the {cc)^(^^'^^Lj) object into 
a physical J/ip by absorbing/emitting soft gluons at the long distance scale. The short 
distance process can be calculated using PQCD in powers of a^, whereas the long distance 
part is treated as a new parameter (OlOg (^'^^^Lj)|0) which characterizes the probability that 
the color-octet {cc)(^^^^Lj) state evolves into a physical J /'ip. By fitting the J /ip production 
at the Tevatron using the usual color-singlet production and the cascades from Xc(l-P) and 
the color-octet contribution, the authors of Ref. |^ determined 

(0|O^(=^5i)|0) = (6.6 ± 2.1) X 10"^ GeV^ (1.3) 

M^+<i^?^ = (2.2±0.5)xlO-GeV' (1.4) 

for Mc = 1.48 GeV. Although the numerical values of two matrix elements {0\Og{^Po)\0) 
and {0\Og{^So)\0) are not separately known in Eq. (1.4), one can still extract some useful 
information from it. Since both of the color octet matrix elements in Eq. (1.4) are positive 
definite, one has 

< (0|C»|'(^^o)|0) < (6.6 ± 1.5) X 10-2 GeV^ (1.5) 
< (^l^t^PoM ^ ^2.2 ± 0.5) X 10-2 GeV\ (1.6) 

These inequalities can provide us with some predictions on various quantities related with 
inclusive J/ ip productions in other high energy processes, which enables us to test the idea 
of color-octet mechanism. 
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Since the color-octet mechanism in heavy quarkonium production is a new idea proposed 
in order to resolve the ip anomaly at the Tevatron, it is important to test this idea in 
other high energy processes with inclusive heavy quarkonium productions. Up to now, the 
following processes have been considered : J/ip production at the Tevatron and fixed target 
experiments f^]) spin alignment of the color-octet produced J/?/) 0, the polar angle 
distribution of the J/ip in the e^e~ annihilations into J/ip + X 0, inclusive J/ip production 
in B meson decays |TD[ and the decays at LEP |]TT|. These processes also depend on 
the aforementioned three color-octet matrix elements in different combinations from (1.4). 
Thus, one can check if the color-octet mechanism provides reasonable agreements between 
PQCD and the experimental data on inclusive J /if) production rates from these processes. 

In the above list of various inclusive J/ip productions at high energies, the J/ip pho- 
toproduction is missing, however. It is the purpose of this work to study the color-octet 
mechanism in the J/ ip photoproduction 0. 

In Sec. we demonstrate how to get the inclusive production rate of a heavy quarkonium 



in the NRQCD formalism of Bodwin, Braaten and Lepage [|1J] . Then, we review briefly the 
7(7 fusion into J/tp + g in the color-singlet model in Sec. p.11 A| . Then in Sec. p.llB| , we consider 
the color-octet subprocesses 

7 + ^7-(cc)C5f or^^pjL^o,^), (1.7) 

which have not been included in previous studies. The size of these color-octet contributions 
to the J/ ip photoproductions are suppressed by relative to the color-singlet contributions, 
but of lower order in a^. This subprocess contributes to the J/ip photoproduction in the 
forward scattering (the elastic peak) with z ^ 1 and ^ 0. These color-octet 2 1 
subprocesses can also contribute to the 2 — > 2 subprocesses through 

^ + g-.(cc)Csl,'^ or ■'Pf^)+g, (1.8) 
7 + (cc)(^5i'^ or ) + g. (1.9) 

These are also resolved photon processes at lower order [0(aa^)] than the color-singlet model 
[0(aa^)] in the perturbation expansion in : although the color-octet contributions are 
suppressed by compared to the color-singlet resolved photon process. Thus, the color- 
octet ^5*0 and ^Pj states can contribute to the elastic peak of the J/V'— photoproduction 
as well as contribute to the resolved photon process. It is quite important to estimate the 
latter and compare with the resolved photon process in the color-singlet model, since it 
is a common statement that J/'?/'— photoproduction is a good place to measure the gluon 
distribution function in a proton. We flnd that the quark contributions are small compared 
to the gluon contribution even if we include (1.9). When one considers Eqs. (1.8) and (1.9), 
one has to consider 

j + g^(cc)C'^'Lf)+g, (1.10) 



^ While we were finishing this work, we received two preprints which discussed the same topic |12| 



|13|. We find our results agree with theirs in the case the direct comparison is possible. 
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although it is expected to be suppressed relative to the usual '-yg fusion color-singlet diagram 
by v^. It is gauge invariant by itself, and thus can be safely neglected if we wish. We 
keep it however, in order to be consistent in the expansion, and make it sure the 
scaling rule works in this case. All of these color-octet 2 — > 2 subprocesses are discussed 
in Sec. [Ill C| . Numerical analyses relevant to the fixed target experiments and HERA are 
performed in Sec. |IV A| . We show that the relations (1.5) and (1.6) yields too large a cross 
section for the J / ijj photoproduction in the forward direction. They also leads to too rapidly 
growing da / dz distribution for high z region compared to the experimental observations. In 
Sec. LV B| , we briefly digress on the B — »• J/ijj + X using the factorization formula derived in 



Ref. [|T0|, and find again that the relations (1.5) and (1.6) overestimate the branching ratio 
for B — s> J/ip + X. All of these seem to indicate that the relations (1.3) and (1.4), especially 
the latter, are probably overestimated by an order of magnitude. This is not surprising at 
all, since the analyses in Ref. employed the leading order calculations for the color-singlet 
parton subprocess for the J/ip hadroproduction. We summarize our results and speculate 
the origins of these overestimates of J/ip photoproductions and B meson decays in Sec. 0. 



II. NRQCD FORMALISM FOR HEAVY QUARKONIUM PRODUCTIONS 

To begin, we consider general methods to get the NRQCD cross section of the process 
a + b —>■ QQC^^'^^L^j'^"'^) H ,where H is the final quarkonium state and QQ{^^~^^L^j'^"'^) 
is the intermediate QQ pair which fragments into a specific heavy quarkonium state in the 
long distance scale. If the on-shell scattering amplitude of the process A{a + b ^ Q + Q) is 
given, we can expand the amplitude in terms of relative momentum q of the quarks inside 
the bound state because the quarks which make the bound state are heavy. Scattering 
amplitude of the process a + b ^ QQ{^^~^^L^j'^"'^){P) is given by 



LzSz siS2 ij 



X (si ; S2 1 SSz) {LLz; SSz \ JJz) (3z3j 1 1 , 8a) 

xA{a + b ^ + q) + Q^{^ - g)), (2.1) 

where the superscript (1,8a) represents the color-singlet or the color-octet configuration of 
the QQ state. After integrating over the relative momentum q, we get 

A{a + b^ QQC'^'L^j''''^)) = ^JcIm', {a + b^ QQC'^'L^j''''^)) , (2.2) 

where 

M'L{a + b^QQC'^'L'^j''''^)) = Yl T.T.i^i-^^2\SSz){LLz;SSz\JJz){3i3j\l,8a) 

LzSz «i«2 ij 

( A{a + b^Q' + Q%=o iL = S), 

x\ e*,iLz)A''ia + b^Q' + Q%=o iL = P),, (2.3) 
[ e;^(Lz)^"^(a + b^Q^ + Q%=o (L = D), 

= y X \q\^^ and C^' = ,7,, x \q\^^+\ (2.4) 

2g0(27r)5 ^ 2q^{2ny ^ ' 
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Here, Cp' is the factor for the relativistic correction of the L state. ^Prom now on, amphtudes 
with Lorentz indices mean 



A^{P,q) 



d 



A{P,q) and A"^(P,g) 



A{P,q). 



(2.5) 



The color SU(3) coefficients are given by 

(3i;3j|l) = Sij/^c and (3i;3j|8a) = ^27;^.. (2.6) 

For the case of color singlet state, we can relate the coefficients to the radial wavefunction 
of the bound state as 



Cs^^\Rsm\ Cp^^lR'Mf and Cn^^Km'- 
Some identities of color matrix trace are useful in squaring the matrix elements. 
Tr(l) = +iV„ 



(2.7) 



Tr(r«r«) = +^"^ 



rj\j.^rparpbrj-iarj-ib^ _ 

Tr(r"r'')Tr(r"r'') = + 



2 ' 

JVi?-l 



Tr (2^oy by fej^c^^a ^ 



iV2-l 



4 ' 



Tr(T'"T'''T'^)Tr(T'T''T''^) = 
Tr(T'"T'''T'^)Tr(T''^T'^T''') = 



(2.8) 



(jV2_l)(Jv2_2) 

4W^ 



At this stage we can derive the explicit form of the matrix element Ai'j^. In general, the 
on-shell amplitude can be expressed as 

(3i3j|l, 8a)A{a + b ^ + q) + Q^{^ - q)) = ui^ + q; Si)0{P, q)v{^ - q; S2), (2.9) 

where O is the matrix relevant to the on shell amplitude. If we introduce the spin projection 
operator Pss^, {P, q) as 



P _ P 

Pss.iP, q)ij = ^{si, S2\SSz)vi{— - q; S2)%(— + q; Si), 

S1S2 ^ ^ 

we can simplify the form of the matrix element Al^ as 



A1'5 = Tr[O(P,0)P5S.(P,0)], 



O'^Pss, + OP 



ss. 



M'n = E ^*apiLz){LW, SS,\JJ,)Tt [O'^^Pss. + O^P^s. + O^P^s. + OP, 

LzSz 



al3 
SSz 



q=0 



(2.10) 



(2.11) 
(2.12) 

(2.13) 



Note that O includes the color coefficient (3i3j|l, 8a) and Pss^ possesses the spin coefficient 
(si; S2\SSz)- Expanding Pss^iP, q) to the second order of the relative momentum q, we get 
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Ao(P,0) 
^o"o(^,0) 

^5.(^,0) 



2V2 



2V2Mc 
1 

2V2Mc 



75(/^ + 2Mq) 



2^2 
1 



1 



/*(5,)(/- + 2Mq), 

[/*(5.)(/P + 2Mq)7" + 7" /*(5.)(/P + 2Mq)] 
1 



-(/^ + 2MQ)(e*"(5,)7'5 + e*^(^,)7 



(2.14) 
(2.15) 
(2.16) 
(2.17) 
(2.18) 

(2.19) 



When L = P, we need further relations to get the correct polarization state of the interme- 
diate state, 



LzSz 



e*^(5,)(lL,;15,|J = J, = 0) 



1 

73 



papP 
M2 



E e*"^''^^e*''(5.)(lL,; 15,| J = 2 J,) = e*"'3(J.), 



P.el(J, 



(2.20) 



where the polarization vector and symmetric polarization tensor have the properties 

P„e-(J,)=0, P„e-^(J,) = 0, e^(J.)=0, e-^(J,) = e^-(J,). (2.21) 
Once the cross section of the on-shell parton level process is calculated, one can expand 



it in factorized forms following BBL [14 



as 



a (a + b ^ {QQ)n 
a[a + b-^ {QQ)n ^H + X 



m: 



Q 



2J+1 ' 
2J+1 ■ 



(2.22) 
(2.23) 



We use a instead of a, as a subprocess cross section, since we will consider 7^ collision where 
the particle b is treated as a parton inside a proton. The index n denotes the intermediate 



QQ state 



25+1 r(l-^ 



which may differ from that of H. The factor multiplied to the H 



production cross section differs from that of QQ state by unity in mass dimension. This 
makes the long range factor coincide with the conventionally normalized wave function of 
the bound state for the color singlet case. We extracted the factor 1/ (2J+ 1) in advance to 
avoid the unnecessary factor after imposing the heavy quark spin symmetry property as 



{0|C>^(3gf'''^)|0) 
3 

(oic'^(3py''^)io) 

2J + 1 



(0|0'^('Po^'''^)|0). 



(2.24) 
(2.25) 
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The 4-fermion operator On with the dimensiond„ are defined as 
dn = 6 rf„. = 8 



OiCPi)=V^(-|L>)x-xt(-|m 



(2.26) 



and dimension-8 operators related to the relativistic correction are 



V'^^X-xM-| L')'^ + h.e 
V'V^xx^fj^X-l)^ B B ^'V + h.c. 



(2.27) 



where 



X^ Dil^ = x\DiIj) - (L'x)V, 



(2.28) 
(2.29) 



and D is the covariant derivative. There are Pauh spinor fields in the previous equations. 
-0 annihilates a heavy quark Q and x creates a heavy antiquark Q. Color and spin indices 
on the fields x have been suppressed. 

Vacuum expectation values of the production operators O^^ and are 



(0|Of |0) = (OlxXV' [Y,Y.\H + X){H + X\\ V^XxlO), 



(2.30) 



(oio^^'^io) = (0|xX^ ( E \QQ{''^'L^}''')){QQ{''^'Lr^)\ ) V^^Olo) 



r(l,8)N 



(2.31) 



= (2J + i)(o|x^/c„v|gQ(^^'+^4«)))(gg(^^'+^L?'«))IV'^/c;x|o) 

= (2J+1)(O|0^'5(25'+1^(1,8))|Q^_ 

The factors /C„ and /C^ are products of a color matrix, a spin matrix and a polynomial in 
the covariant derivative D and other fields , which are same with those of On- According to 
the heavy quark spin symmetry, 2S+i^0^>8) state has the same properties with another state 
2«'+iL?'^Vith the same L) except for the mj multiplicity factor 2 J + 1, which appears in 
the last equation. 

Let us explain the process to derive the short distance coefficients Regardless of the 
onshell scattering amplitude a + h ^ Q -\- Q, the intermediate bound state QQ{^^^^ L^^'^^) 
production cross sections have common factor (0|C»^'5|o). We can obtain the matrix elements 
by using the intermediate state ket 



\QQ{ 



25+1 r(l,8a) 



){P)) 



E EE 

LzSz «i«2 ij 



(27r)32g0 



2M, 



Q 



LL 



X (si;s2|-55z)(LLz;-S5z|JJz)(3z3j|l,8a)|gX^ + ?)^'(^-?)) (2-32) 



as 

2J + 1 

where 



Cl X Cn, (2.33) 



_ J 2Nc (color singlet) , . 

- I iV2 _ 1 (color octet) " ^^"^^^ 

The cross section of the process a{pi) + b{p2) — > {QQ)n{P) (with n representing the partial 
wave and the color quantum numbers of QQ) is given by 

a{a{p,)+b{p2)^{QQUP)) 

^ a'^ X Cl X Cn 

where 

K = 7^?<^(^ - M|)EI-Ml(« + b - (gg)n)r, (2.36) 

n represents the partial wave (^^~^^Lj) and the color quantum numbers of QQ and s is the 
invariant mass of the initial particles a and b. Then, the bound state cross section and the 
short distance coefficients F„ are given by 

a{aip,) + b{p,) {QQ)n ^H + X) = — x = x (2.37) 

F„ = a; X M^Q--\ (2.38) 

In case of 7p scattering, we should convolute the above result with the parton structure 
functions to get the cross section : 

^(a(pi) + 6(P2) - {QQ)n ^H + X) = ^Y. <ib) X ^^jff' (2-39) 

where 

f^'nib) ^ j dxfb/p{x)a'n 



n 



16CM 



^/Vp(^1_,,..,.EI-^l(« + ^^ (W)n)r. (2.40) 



x=4M^/ 



For the case of 2 — > 2 process, we need to modify the formulae only a little. The quarkonium 
H photoproduction cross section via 2 — > 2 subprocess a{pi)+b{P2)/p — > {QQ)n{P)+c{P3) 
H + X, where 6 is a parton of the initial proton is given by 
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da{a{p,) + b{P2)/p ^ (QQUP) + c{ps) ^ H + X) 



■dzdP^. (2.41) 



For example, if we consider the J/ip production via ^5*0*^^ '^Sf'\ ^Pq'\ ^Pi^^ and ^^2*^ 
intermediate states, we get 

V3p(8),^ (0|O^(^Pf)|0)^ 



Mr 



j:^'cpj 

J 



2J+ 1 



{0\O^CS^''^)\0) X (a'(iS^'^) + ^'(^5}'^)) + (0|C'^(3Po^'^)|0) X ^a'(3pj^)^ 

J 



(2.42) 



III. J/i; PHOTOPRODUCTION SUBPROCESSES 

A. Color-singlet contributions 

The inelastic J/?/'— photoproduction has long been studied in the framework of PQCD 
and the color-singlet model ||15|. The lowest order subprocess at the parton level for 



"f + p^J/ip + Xis the 7— gluon fusion at the short distance scale (Fig. 1), 

j + g^(cc)(^S['^)+g, (3.1) 

followed by the long distance process 

(cc)(35«)^ JM (3.2) 

at the order of 0{aa'^v^) in the nonrelativistic limit. Thus, the production cross sec- 
tion is proportional to the gluon distribution inside the proton. This is why the 
J/'?/'— photoproduction has been advocated as a clean probe for the gluon structure function 
of a proton in the color-singlet model. Without further details, we show the lowest order 
color-singlet contribution to J/tp photoproduction through 7— gluon fusion in the nonrela- 
tivistic limit : 



s\s - AM^y + t\t - AM^y + u{u - 4M, 



where 



z 
s 



St- I lab pj^.fc) 



{k + q,y = xs, (3.4) 



t = {P-ky = {z-l)s. 
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The overall normalization TVi is defined as 



M = - (47r«.)'(47r«)e,^ MlGi{J/ij). (3.5) 

y 

The parameter Gi{J/ip), which is defined as 

Gi{J/ij) = (3.6) 

in the NRQCD, is proportional to the probability that a color-singlet cc pair in the 'S['^ 
partial wave state to form a physical J/ip state. It is related with the leptonic decay via 

r(J/^ l+r) = - nela"^ Gi{J/^), (3.7) 
3 

to the lowest order in a^. From the measured leptonic decay rate of J/ip, one can extract 

Gi(J/^) ^ 106 MeV, (3.8) 

Including the radiative corrections of 0{as) with as{Mc) = 0.27, it is increased to ~ 184 
MeV. Relativistic corrections tend to increase Gi{J/ip) further to ~ 195 MeV |10[ . 
The partonic cross section for 7 + a— »J/-?/^ + 6is given by 

§ = TTT^ El-^(7 + a - -//^ + (3.9) 
dt Ib'KS'^ ^ 

The double differential cross section is 

5^(0^ + ^ JmP,enX) ^ ^ ^ T.\M\XsA (3.10) 

where 

s 1 



s zs 



p2 

Ml+ 



^ ■ l-z 



(3.11) 



One has the following constraints for x, z, t and 



2 . 



^<x<l, 3.12 
s 

-(i-Ml)<i(=t)<0, (3.13) 

The z and distributions can be obtained in the following manner : 

dP|, (3.15) 



dz Jo dzdPj' 
da /-^max d^a 



Zrr 



z„ 



\;[s + Ml + ^{s-MlY-AsP(), (3.17) 



2s 
1 

2i 
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There are two kinds of corrections to the lowest order result in the color-singlet model, 
(3.1) : the relativistic corrections of O(u^) and the PQCD radiative corrections of 0{as) 
relative to the lowest order result shown in (3.1). We briefly summarize both types of 
corrections in the rest of this subsection, since they have to be included in principle for 
consistency, when one includes the color-octet mechanism in many cases. 



The relativistic corrections to the 7— gluon fusion was studied by Jung et al. []T5 
They found that relativistic corrections are important for high z > 0.9 at EMC energy 
i^/^ip — ^^-"^ GeV). Since it mainly affects the high z region only, we neglect the relativistic 
corrections, keeping in mind that it enhances the cross section at large z > 0.9. 

The radiative corrections to the J/tp photoproduction is rather important in practice. 
This calculation has been done recently in Ref. [|l^], and the scale dependence of the lowest 
order result {Q^ in the structure function in Eq. (3.10) ) becomes considerably reduced. 
For EMC energy region, the K factor is rather large, K ^ 2. For HERA, it depends on 
the cuts in z and P^. We include the radiative corrections in terms of a K factor suitable 
to the energy range we consider. Another consequence of the radiative corrections to the 
color-singlet J/ip photoproduction is that the PQCD becomes out of control for z > 0.9 
at EMC energy. For HERA, one gets reasonable results in PQCD when one imposes the 
following cuts in z and : z < 0.8 and P^ > 1 GeV^ . Thus, it is not too much of sense to 
talk about the z or px distributions for such z region in PQCD. One has to introduce cuts 



in z as well as in px- Following the Ref. |jT6[, we adopt the following sets of cuts 



z < 0.9, for EMC, (3.19) 
0.2<z < 0.8 for HERA. (3.20) 

At HERA energies, the lower cut in z{z > 0.2) is employed in order to reduce backgrounds 
from the resolved photon process and the b decays into J/V'. For these cuts, the K factor is 
approximately ~ 1.8 both at HERA and the fixed target experiments. We include these 
radiative corrections to the subprocess (3.1) in Sec. [IV A| by setting K ~ 1.8. 



B. Color-octet contributions to 2 ^ 1 subprocesses 



Let us consider color-octet contributions to the 2^1 subprocesses via 



(3.21) 



followed by {cc)s fragmenting into J/tp with emission of soft gluons. This subprocess occurs 
at 0(aasf''). Here, a,b are color indices for the initial gluon and the final color-octet cc 
state, and we are interested in5' = L = J = and S' = L = 1, J = 0, 1, 2. There are 2 
diagrams representing the vertex, as shown in Fig. 2. Here we consider the process where 
only the gluon is off-shell. Following previous convention, we first write the matrix O related 
to this effective vertices. 



0{P, q) 



72 



ab 



I + g - A;)2 - Ml 



f + ^- ^ + Me 



(3.22) 
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With this matrix O we can derive the effective vertices for the jglcc 



,2S+lr(8) 



Ly' as 



7W'(35f )) = 0, 



V2eecgsS'''' 



g^ - 4M2 i 



(3.23) 
(3.24) 
(3.25) 



M2(c/2 _4M2) 
X (^gh""''^ + 2k^ 
IQee^gsd''^ 



^ ^ ' eaWeMPr, (3.26) 



ig' - 4M2) 



2k 



g^ - ml J 

^2 - 4M2 



e«^(J,)e;ie^ (3.27) 



Since J/ip can be produced via the 2—^1 subprocesses with this effective vertices we can 
obtain the 2^1 color octet contribution by using the following average squared amplitudes 
asQ 



3p(8 



Y.\M'{'P? 
Y.\M'{'P^' 



2{eecgs 

0, 
6 . 



M2 



M2 



[eccgs 



[eccgs 



(3.28) 
(3.29) 

(3.30) 

(3.31) 
(3.32) 



The photoproduction cross section of the J/ ip production via 2—^1 process can be obtained 
from (|2.42|), assuming heavy quark spin symmetry : 



a 



(7 + 



P 



fcc)(«) 



77r(eecgf^)^ 
64M5 



I {0\O^{^Pi^^)\0) (0|C>^(^g^'^)|0) ' 



(3.33) 



Since a oc 5(1 — 2;), this 2^1 color-octet subprocesses contribute to the elastic peak in the 
J/'?/'— photoproductions. It is timely to recall that the color-singlet model with relativistic 



corrections still underestimate the cross section for z > 0.9 by an appreciable amount [15 



As 2 — i> 1, the final state gluon in the 7— gluon fusion becomes softer and softer, although 
this does not cause any infrared divergence in the transition matrix element. Therefore, 



20ur results agree with those obtained in Refs. 
the invariant matrix is different from theirs. 



. Note, however, that our convention of 
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it would be more meaningful to factorize the effect of this final soft gluon into the color- 
octet matrix elements, {Og{^So)) and {Og(^Pj)). The color-octet ^5*0 and ^Pj states might 
reduce the gap between the color-singlet prediction and the experimental value of da/dz for 
0.9 < ;z < 1. 



C. Color-octet contributions to the 2^2 subprocesses 

The color-octet 2 — 1 subprocess (3.21) considered in the previous subsection not only 
contribute to the elastic peak of the J/ip photoproduction, but it also contributes to the 
resolved photon processes at O(aa^v^), shown in Fig. 3, where the initial partons can be 
either gluon or light quarks (g = u, d, s). These diagrams are suppressed by v"^ but enhanced 
by 1/as, relative to the resolved photon process in the color-singlet model. Also, the jet 
structures from Fig. 3 are different from that from the resolved photon process in the color- 
singlet model. They can enhance the high-pj- J/^'s, which might be relevant to HERA 
energy. This can be a background to the determination of gluon distribution function of a 
proton, if the cross section is appreciable. The resolved photon process in the color-singlet 
model is dominant over the 7— gluon fusion in the lower z region, z < 0.2, and it can be 
discarded by a suitable cut on z. Since the color-octet contribution to the resolved photon 
process has not been studied in the literature, we address this issue in the following. When 
one considers Fig. 3, one has to include Fig. 1 with {cc)s simultaneously, since both are the 
same order of 0{aa1v'^). This diagram is the same as the color-singlet case except for the 
color factor of the (cc) state. 

It is straightforward, although lengthy, to calculate the amplitudes for the processes 
shown in Figs. 1 and 3. Using REDUCE in order to the spinor algebra in a symbolic 
manner, we can get the averaged Ai squared for various 2^2 processes. Since the full 
expressions are rather involved, they are shown in Appendix A. 

Another color-octet {cc){^S[^^) contribution to the J/-?/^— photoproductions comes from 
the Compton scattering type subprocesses (see Fig. 4) : 

l{k,e)+q{p,) ^ {cc)est\P,e*) + q{p,), (3.34) 

where P and e* are the four momentum and the polarization vector of the ^5*1 color-octet 
state, and a is its color index. This subprocess, if important, can be a background to the 
determination of the gluon distribution function in a proton, since it is initiated by light 
quarks. /^From the naive power counting, however, we infer this subprocess occurs at 0{aal) 
in the coupling constant expansion, and also suppressed by compared to the color-singlet 
contribution (3.1) due to its color-octet nature. Thus, this subprocess is expected to be 
negligible. 

One can actually quantify this argument by explicitly evaluating the Feynman diagrams 
shown in Fig. 4. The effective vertex for qq — * {cc){^Si^^^) is given by (Fig. 5) Q 

M'iqipMP2) - (cc)('^r)) = ^ vip2)rT'^uip,) e;(pi +p2,^.), (3.35) 

where e* is the polarization of the produced spin-1 color octet object. Using this effective 
vertex, one can calculate the amplitude for the Feynman diagrams shown in Fig. 4 : 
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A 



(A; + pi)2 -M2 
(pi - P + M, 



{pi - P)2 - M2 



u{p{). (3.36) 



where eCq is the electric charge of the hght quark inside proton(g = u,d,s). The averaged 
Ai"^ for the color-octet ^Si state is given by 



El-M'(7g-(cc)(='5f))g)| 



u 

+ - + 



3M2 " " ~M s su 



(3.37) 



This completes our discussions on the color-octet 2^2 subprocess for J/ip photoproduc- 
tions. 



IV. NUMERICAL RESULTS 

A. J/?/; photoproductions at fixed targets and HERA 

Now, we are ready to show the numerical results using the analytic expressions obtained 
in the previous section. Let us first summarize the input parameters and the structure 
functions we will use in the following. The results are quite sensitive to the numerical values 
of Us and ttlc and the factorization scale Q. We shall use ^^(M^) = 0.3, ttlc = 1.48 GeV 
and = (2mc)^. For the structure functions, we use the most recent ones, MRS A |0 and 
CTEQ3M [0 , which incorporate the new data from HERA , on the lepton asymmetry 



in 1^— boson production [^] and on the difference in Drell-Yan cross sections from proton 



and neutron targets |^T|. For the 2 — >■ 1, we show results using both structure functions. 
For the 2 — 2 case, we show the results with the CTEQ3M structure functions only, since 
the MRSA structure functions yield more or less the same results within ~ 10% or so. 
Let us first consider the J/^p photoproduction via the color-octet 2-^1 subprocess 



considered in Sec. [1IIB| . Since the subprocess cross section (3.33) vanishes except at z = 1, 
one can infer that it contributes to the J/ip photoproductions in the forward direction 
[z ~ 1, ^ 0). In Figs. 6 (a) and (b), we show the J/tp photoproduction cross section in the 
forward direction as well as the data from the fixed target experiments and the preliminary 
data from HI at HERA, respectively. In each case, the upper and the lower curves define the 
region allowed by the relation (1.4) for two color-octet matrix elements, {O\0g Si)\0) and 
(OjOg (^Po)|0). In case of fixed target experiments, it is usually characterized by 2 > 0.9, 
with the remainder being associated with the inelastic J/ip photoproduction. According 
to this criterion, the experimental value of crexp(7 + p ^ J/ip + X) contains contributions 
from inelastic production of J/ip^s. Thus, the data should lie above the predictions from the 
color-octet 2 — > 1 subprocess, (3.21). Fig. 6 (a) shows that the situation is opposite to this 
expectation. Color-octet contributions are larger than the data, which indicates that the 
numerical values of the color-octet matrix elements are probably too large. At HERA, one 
has the elastic J/ip photoproduction data, which can be identified with the color-octet 2^1 
subprocess. By saturating the relation relation (1.4) by either color-octet matrix element, we 
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get the J /if) photoproduction cross section in the forward direction (Fig. 6 (b)). We observe 
again that the color-octet contribution with (1.4) overestimates the cross section by a large 
amount. This disagreement can arise from two sources : (i) the radiative corrections to 

— >■ J /ip + X, which were ignored in Ref. |^ is important, and/or (ii) the heavy quark spin 
symmetry for (OlOg (^Pj)|0) ~ (2J + 1) (OlOg (^Po)|0) may not be a good approximation. 
Although the heavy quark spin symmetry relation is used quite often in heavy quarkonium 
physics, it may be violated by a considerable amount |T^. 

Recently, Fleming et al. performed the fit to the available fixed target experiments 
and the HERA data independently, and found that 

(0|Og^(i5o)|0) + ^ (0|Og^(^Po)|0) = (0.020 ± 0.001) GeV^ (4.1) 

using the MRSA^, and CTEQ3M structure functions with a,(2Mj = 0.26 and = 
1.5 GeV. This determination is not compatible with the relation (1.4), since the resulting 
(Og(^Po)) is negative. This is another way to say that the determination of the color-octet 
matrix elements from the J/tp productions at the Tevatron may not be that reliable. In fact, 
this is not very surprising, since the radiative corrections to the lowest-order color-singlet 
contributions to J/ip hadronic productions are not included yet. 

Next, we consider J/ip photoproductions through 2 — >• 2 parton-level subprocesses. As 
discussed at the end of Sec. [Ill AJ , the PQCD corrections to the lowest order 7-1-5' J /ifj + g 
is not under proper control for z > 0.9. Therefore, we impose a cut z < 0.9 at EMC energy, 
^s^p = 14.7 GeV, and at HERA with y/s^ = 100 GeV, we impose cuts on z and P^ [|16| : 



0.2< z < 0.8, P4 > 1 GeV\ 

In both cases, we set K ~ 1.8. 

In Figs. 7 (a) and (b), we show the da/dz distributions of J/iJj at EMC (NMC) and HERA 
along with the corresponding data. In both cases, the color-octet ^5*1 contribution (Compton 
scattering type) is negligible in most regions of z, and thus can be safely neglected. The 
thick dashed and the thin dashed curves correspond to the cases where the relation (1.4) 
is saturated by Og(^Pj) and Og(^5'o), respectively. The thick and the thin solid curves 
represent the sum of the color-singlet and the color-octet contributions, in case that the 
relation (1.4) is saturated by (0|Og (^S'o)|0) and (0|Og (^Po)|0), respectively. In either case, 
we observe that the color-octet ^Sq and ^Pj contributions begin to dominate the color-singlet 
contributions for z > 0.6, and become too large for high z region considering we have not 
added the enhancements at high z due to the relativistic corrections. Thus, this behavior 
of rapid growing at high z does not agree with the data points at EMC and HERA, if we 
adopt the determination (1.4) by Cho and Leibovich [^]. 

In Figs. 8 (a) and (b), we show the P^ distributions of J ftp at EMC and HERA, 
respectively. We find that the color-octet contributions through 2—^2 subprocesses become 
dominant over the color-singlet contributions for most P^ region. Also, the color-octet 
contributions from ^Sq and ^Pj are more important than the charm quark fragmentation 



considered by Godbole et al. |Q. However, this situation may be due to the overestimated 
color-octet matrix elements as alluded in the previous paragraph. 

In Fig. 9, we show the inelastic J/ip photoproduction cross section as a function of y/s^ 



with the cut, z < 0.8 and P^ > 1 GeV^. Again, the color-octet ^5*1 contribution is too small. 
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and thus not shown in the Figure. Here again, the color-octet ^5*0 and ^Pj contributions via 
2 — i> 2 subprocesses (Fig. 3) dominate the color-singlet contribution, if the relation (1.4) is 
imposed. Although the total seems to be in reasonable agreement with the preliminary HI 
data, direct comparison may be meaningful only if the cascade J/ip^s from b decays have 
been subtracted out. There are also considerable amount of uncertainties coming from Mc 
and as- Therefore, it is sufficient to say that the color-octet ^5*0 and ^Pj state dominate the 
singlet contribution to J/ip photoproduction, if the relation (1.4) is imposed. 



B. Digression on S — > J/tp + X 

Finding that the color-octet matrix elements determined from J/tp productions at the 
Tevatron seem to be too large in case of J/ ip photoproductions at fixed target experiments 
and HERA, it is timely to consider the color-octet contributions to inclusive B meson decays 
into J/ip + X again. In Ref . |TD[ , a new factorization formula is derived foT B —>■ J/ip + X : 



{o\oi'*('p,m ,^ . ^ /, . »Mi 



with 



Gl\ I. ml 

2 



Two numbers C+(Mb) ^ 0.87 and C_(Mfe) ^ 1.34 are the Wilson coefficients of the |A5| = 1 
effective weak Hamiltonian. Using the relation (1.4), we estimate the above branching ratio 
to be (for ^.(M^) = 0.28 in Ref. |ig) 



(0.42% X 12.8) < B[B J/iP + X) < (0.42% x 13.8) (4.4) 
which is larger than the recent CLEO data by an order of magnitude 

5cxp(5 J/ip + X) = (0.80 ± 0.08)%. (4.6) 



^Even if we use the new determination (4.1) by Fleming et al. we still get a large branching 
ratio : 

(0.42% X 3.45) < B{B J/iP + X) < (0.42% x 5.45), (4.5) 
although the discrepancy gets milder than the case (4.4). 
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The situation is the same for B ip' + X. This is problematic, unless this large color- 
octet contributions are canceled by the color-singlet contributions of higher order in 0{as) 
which were not included in Ref. [|l^]. If there are no such fortuitous cancelations among 
various color-octet and the color-singlet contributions, this disaster could be attributed to 
the relation (1.4) being too large compared to the naive velocity scaling rule in NRQCD, as 
noticed in Ref. [Q. It seems to be crucial to include the higher order corrections of 0{a'^) 
for the color-singlet J/ip productions at the Tevatron, which is still lacking in the literature. 



V. CONCLUSION 

In summary, we considered the color-octet contributions to (i) the subprocess 7 -|- p — »• 
J/ip + X through 7(yf — {cc)%{^Sq and ^Pj) and the subsequent evolution of (cc)8 into a 
physical J /ip with ^ 1 and ~ 0, (ii) the subprocesses 7 + (y'(or q) {cc)%{^Sq or ^Pj) + 
g{oi q). These are compared with (i) the measured J/ip photoproduction cross section in 
the forward direction, and (ii) the z distributions of J/ip at EMC and HERA, and the 
preliminary result on the inelastic J/ip photoproduction total cross section at HERA. One 
finds that the relation (1.4) color-octet lead to too large contributions of the color-octet ^Sq 
and ^Pj states to the above observables. Especially, the first two observables contradict the 
observation. This is also against the naive expectation that the color-octet contribution may 
not be prominent as in the case of the J/ hadroproductions, since they are suppressed by 

(although enhanced by one power of a^) relative to the color-singlet contribution. It is 
also pointed out that the same is true of the process B ^ J/ip + X, in which the relation 
(1.4) predicts its branching ratio to be too large by an order of magnitude compared with 
the data. 

Therefore, one may conclude that the color-octet matrix elements involving ccs{^So, ^Pj) 
might be overestimated by an order of magnitude. Since the relation (1.4) has been extracted 
by fitting the J/ ip production at the Tevatron to the lowest order color-singlet and the color- 
octet contributions, it may be changed when one considers the radiative corrections to the 
lowest order color-singlet contributions. 

While we were finishing our work, there appeared two papers considering the same subject 



|T^. Our results agree with these two works where they overlap. 
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APPENDIX A: THE INVARIANT AMPLITUDE SQUARED FOR THE 2^2 
SUBPROCESSES CONSIDERED IN SEC. III.C. 

In this Appendix, we give explicit expressions for the invariant amphtude squared for the 
color-octet 2 — > 2 subprocesses shown in Figs. 1 and 3. The results were obtained using the 
symbolic manipulations with the aid of REDUCE package. We checked that the following 
formulae in two independent ways : the helicity amplitude and the covariant density matrix 
methods in order to do the summations over gluon polarizations. Now both methods yield 
the same results listed below : 



cc 



,25+1 r (8) 



167rs2 (2J+l)Me 



3su 



3o(8)^ 



2t(s + t)2(t + M)2(M + s)2 

su 



(Al) 



+ + u^ + (2M, 

s'^{i + uf [2iu{i + u) - 3s2(2Me)^ + isuf + 9(2Mc)*(t + sf{s + uf{u + i) 
W{i + sf [2is{i +s)- 3n^(2Mc)^ + isuf 
+i^{s + uf + s){i + u) + 2i{s + u) + 2{f + u^)'] 

I 2Mli{s + t)\i + u)\u + sY 

3 {2McfP |s^(t + uf (u + {2Mcfy + u^it + sf (s + (2M, 
+ (2M,)2sV |(£ {I -s)u)\ {i +sf{{t + uf + {i 

+ (2Mj2(s2 _ ^2^12 ^ u){t + s) - 2suf 
+2isu ((2Mc)^ - suY (s^(t + uf + u^{i + sf + P{s + uf) 



(A2) 



— UjS 



I 2Ml{s\if{Uu)\u\sf 



(A3) 
(A4) 



Mli{s^if{i+uf{u + sf 



t^(s + M){2(s2 + u2) + 115^1 



^ The original expressions (Al)-(A4) are modified, because of maltreatment of gluon polarizations 
of Fig. 3 with three gluon vertices. The numerical results based on the correct expressions shown 
here were already given in the Ref. |23|, and so we do not reproduce them here. 

JL would like to thank Avto Tkabladze for pointing out this error and checking that the above 
formulae are correct. JL also would like to thank B. Kniehl and G. Kramer for using the corrected 
formulae for their analysis in Ref. |23| where they reproduced exactly the same numerical results 
as in Ref p|. 
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/|2 



+i\s + u) [2{s* + ii^) + 61sm(s2 + M^) + 75s^u^} 
+?s?'},(s + ?-},) {25(s^ + ii^) + 63su{s^ + m^) + 93sV} 

+7s'^u^(s + u)(s'^ + SU + u^Y 

(7 + ? - (cc) (^4')) + ,) = -^(ee.y,^)^i^, 



(7 + ? 



El-^f (7 + ?- (cc) ('Pf ) + g) 



32/ 2\2 



i{s + m)4M2 ' 



(A5) 
(A6) 
(A7) 

(A8) 



Here, the Cc = 2/3 and we have summed over the electric charges of light quarks (g = u, d, s) 
in the above expressions, assuming rUg = 0. 
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FIGURE CAPTIONS 



Fig.l Feynman diagrams for the color-singlet and the color-octet subprocess 7 + fif — > 
{cc),,s{'S,)+g. 

Fig. 2 Feynman diagrams for the color-octet subprocess J + g ^ {ccjsi^So or ^Pj). 

Fig. 3 Feynman diagrams for the color-octet contribution to the resolved photon 7 -|- 
g(oT q) (cc)8(^So or ^Pj) + g(or q). 

Fig. 4 Feynman diagrams for the color-octet 2^2 subprocess, 7 + Q' — > {cc)s{^Si) + q with 
q = u,d, s. 

Fig. 5 A Feynman diagram for qq — > (cc)s{^Si). 

Fig. 6(a) The cross sections for 7+p^ J/ip + X in the forward direction at the fixed 
target experiments as a function of E^. The solid and the dashed curves were obtained using 
the CTEQ3M and the MRSA structure functions. Here, TOT^ is the ^Sq^^ saturated curve 
and TOTp is the ^P^ saturated one. 

Fig. 6(b) The cross sections for 7 + p J/ip + X in the forward direction at HERA 
as a function of the square root of s^p. The solid and the dashed curves were obtained using 
the CTEQ3M and the MRSA structure functions. Here, TOT^ is the -^Sq^^ saturated curve 
and TOTp is the ^Pf'^ saturated one. 

Fig. 7(a) The differential cross sections da/dz ior ^ + p ^ J/t/j + X at EMC as a function 

of ^ = Eji^/E^. The singlet contributions are in the thick dotted curve, the color-octet ^5*0 
contributions in the thick dashed curve (with (0^(i5o)) = 6.6 X 10-2 GeV'^), and the color- 
octet ^Pj contributions in the thin dashed curve (with {Ot{^Pj))lMl = 2.2 X 10-2 GeV=^). 
The total is shown in the solid curve. The relation (1.4) allows the region between two solid 
curves. Here, TOT^ is the ^S'f' saturated curve and TOTp is the ^Pf' saturated one. 

Fig. 7(b) The differential cross sections da/dz for ^ + p J / ip + X dX HERA as a function 
oi z = Eji,^/E^. The singlet contributions are in the thick dotted curve, the color-octet ^Sq 
contributions in the thick dashed curve (with {Oti^So)) = 6.6 X 10-2 GeV^), and the color- 
octet ^Pj contributions in the thin dashed curve (with {Ot{^Pj))/M^ = 2.2 x lO'^ GeV^). 
The total is shown in the solid curve. The relation (1.4) allows the region between two solid 
curves. Here, TOT^ is the ^5*0^^ saturated curve and TOTp is the ^Pj'^ saturated one. 

Fig. 8(a) The differential cross sections da/dP^ for 7 + p J/ip + X at HERA as a 
function of P^. The singlet contributions in the thick dotted curve, the color-octet con- 
tributions in the thick dashed curve (with {Oti^So)) = 6.6 X 10-2 GeV^), and the color-octet 
^Pj contributions in the thin dashed curve (with {Ot{^Pj))/M^ = 2.2 x 10^2 GeV^). The 
total is shown in the solid curve. The relation (1.4) allows the region between two solid 
curves. Here, TOT^ is the ^5*0^^ saturated curve and TOTp is the ^Pj"^ saturated one. 
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Fig. 8(b) The differential cross sections da/dP^ for 7 + p ^ J/ip + X at HERA as a func- 
tion of of J/ip. The singlet contributions in the thick dotted curve, the color-octet ^Sq 
contributions in the thick dashed curve (with {Oti^So)) = 6.6 X 10-2 GeV^), and the color- 
octet ^Pj contributions in the thin dashed curve (with {Ot{^Pj))/M^ = 2.2 x lO'^ GeV^). 
The total is shown in the solid curve. The relation (1.4) allows the region between two solid 
curves. Here, TOT^ is the ^5*0^^ saturated curve and TOTp is the ^Pjf^ saturated one. 

Fig. 9 Total inelastic J/ip photoproduction cross section for z < 0.8 as a function of the 
square root of Sjp. The singlet contributions in the thick dotted curve, the color-octet ^-S'o 
contributions in the thick dashed curve (with {OtCSo)) = 6.6 X 10-2 GeV=^), and the color- 
octet contributions in the thin dashed curve (with {Ot{^Pj))/M^ = 2.2 x 10"^ QeV^). 
The total is shown in the solid curve. The relation (1.4) allows the region between two solid 
curves. Here, TOT^ is the ^^o^"* saturated curve and TOTp is the ^P^f^ saturated one. 
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